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Abstract 

The charge and spin sectors, which are intimately coupled to the fermionic one, of the t-t' Hubbard model have been 
computed self-consistently within the two-pole approximation. The relevant unknown correlators appearing in the 
causal bosonic propagators have been computed by enforcing the constraints dictated by the hydrodynamics and 
the algebra of the composite operators coming into play. The proposed scheme of approximation extends previous 
calculations made for the fermionic sector of the t-t' Hubbard model and the bosonic sector of the Hubbard model, 
which showed to be very effective to describe the overdoped region of cuprates (the former) and the magnetic 
response of their parent compounds (the latter). 
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Quite recently, we have shown how it is possible 
to capture, by means of the t-t' Hubbard model, the 
single-particle properties of cuprate materials in their 
overdoped region within an approximation scheme 
that has ingredients like: a two-pole reduction of the 
fermionic retarded propagator, the use of composite 
operators and the implementation of algebraic con- 
straints [1,2]. On the other hand, we have also shown 
how it is possible to reproduce the magnetic response 
of the parent compounds of the cuprates within a sim- 
ilar approach applied to the relevant bosonic causal 
Green's functions [3,4]. Both of these schemes have 
proved to be very reliable through many positive 
comparisons to existing numerical data [1,3]. 
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Along this line, we decided to study the possibility to 
extend the latter approach to the t-t' Hubbard model 
in order to provide further elements of analysis to our 
study of the cuprate properties and, in particular, of 
the magnetic properties in the underdoped region. The 
2D t-t' Hubbard model is described by the Hamiltonian 



H = ^(tij-^i j )c t (i,t) C (j,t) + (/^n T (i)n 1 (i) (1) 

■J i 



We use the standard notation: c(i), c'(i) are annihila- 
tion and creation operators of electrons in the spinorial 
notation; i stays for the lattice vector and i = (i,i); 
fi is the chemical potential; £ij denotes the transfer in- 
tegral; U is the screened Coulomb potential; n a (i) = 
cl(i)co-(i) is the charge density of electrons at the site 
i with spin a. For a cubic lattice the hopping matrix 
has the form £ ij = —Ata^ — 4t'/3ij, where a.ij and j3ij 
are the first and second neighbor projection operators, 
respectively, 
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a « = 2¥ E e ik<fli ^ j) [cos(M + cos(fc y )] (2) 



= ^ E e ' k<iii cos (M cos ( k y) 
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We choose as fermionic basis the following doublet 




(4) 





— C \l)<7pC (t) 


+ c (^Jct m c(jJ — AC 


(i)<7,iC a (i) 








(7) 




= ^(^CTpC™ 5 ^ 


)+c^ {i)a IM c(i)-c^{ 


;i)<r M c^(i) 




-C t/3 (i)cT p C a (i) 




(8) 




= C t (i)a M)? "(i) 


- V Hi)a,c a (i)+ V ^(i 






- c ta (i)a M ?7(i) 




(9) 








i<7 (t c(i) 








(10) 



where = [1 — n(i)]c(i) and r)(i) = n(i)c(i) are 
the Hubbard operators, and n(i) — ^2 a n„(i). In the 
two-pole approximation [1] the retarded GF G(i,j) = 
{R[ip(i)ip^ (j)]) satisfies the equation 



In the two-pole approximation the causal GF 
G (M) (M') = (T[N^\i)N^(j)]) satisfies the equa- 
tion 



[w-e(k)]G(fc,w) =/(k) 



(5) 



(k)]G (M) (k,a;) = J (M) (k) (11) 



where /(k) = F.T.({i>(i,t),^Q,t)}) and e(k) = 
F.T.({i^^-,^ (j,t)})r 1 (k)- the symbol F.T. de- 
notes the Fourier transform. In the paramagnetic 
phase the energy matrix e(k) depends on the follow- 
ing set of internal parameters: p, A = — 

(^(0^(0). A' = (€"(t)ft(t)) - faWW). 

p = «(i)n M (i)>/4 - <[ CT (i) Ci (i)]« C {(i) c |(i)), p' = 

(n^(i)n M (i))/4 - ([c T (i)c i (i)] /3 cj(i)c|(i)), which must 
be self-consistently determined. Given an operator 
C(i), we are using the notation C 7 W = Ej 7ijC(j> 
and C 7A (i) = Eji 7ijAjiC(l, t) with 7 , A = a, (3. The op- 
erator n M (i) = c' (z)<7 M c(i) [<r M = (1, a), a are the Pauli 
matrices] is the charge (p = 0) and spin (p = 1,2,3) 
density operator. The local algebra satisfied by the 
fermionic field (4) imposes the constraint (£(i)»7^(i)) = 
0: this equation allows us to solve self-consistently the 
fermionic sector [1]. 

We consider then the composite bosonic field [3] 



N M (i) 




(6) 



where p M (i) = p'^i) + rp£(i), p'^i) = c t (i)a M c a (i) - 
c at (i)cr M c(i) and p'^i) = c t (i)a M c' 3 (i) - c' 31 ' (i)ff M c(i) 
with r = t'/t. The equation of motion of n M (i) reads 
as iJ^n M (i) = — 4ip M (i), whereas that of p^(i) reads as 



-4tlp(i) + Ukfj_(i). According to the fact 



'St 

that r is usually chosen of the order 10 , we have 
decided to neglect terms of the order r 2 and higher, 
as they will practically give no relevant contributions 
to the dynamics. In this case, we have i M (i) = I'^i) + 
2tI'^(i) and K M (i) = + tk'^(i) with 



where I M {k) = F.T.([N M (i,t), N ( ^Q,t)]),m M (k) 



F.T.{[ 



; 8Af(") (i,t) 
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, N (")t (j, i)]> and e (,i) (k) = (k) [7^ (k)] " x . 



As it can be easily verified, in the paramagnetic 
phase the normalization matrix I^' does not depend 
on the index p; charge and spin operators have the 
same weight. The two matrices 7^ and are 
off-diagonal and diagonal, respectively, and have the 
following entries: 



/&> (k) = 4 [1 - a (k)] C& + 4 [1 - /?(k)] tCL (12) 

mtf>(k) = -4*7i afc (k) (13) 

m&>(k) = -4*J,„ P( .(k) + tf/^k) (14) 

C? c = (c^i^i)) (15) 

Wk) = F.T.(MM),pt(j,*)]) (16) 

7 K ^ M (k) = F.T.([K M (i,t),pt(j,t)]) (17) 



The energy matrix (k) has off-diagonal form with 



entries: e^(k) = -At and e£> (k) = m^(k)/7}^(k). 
Quite lengthy calculations shows that the energy spec- 

tra [£<">(k) = (-)" V / -4tm^ ) (k)/7^ ) (k) n = 1,2] 

and the spectral weights (which can be expressed in 
terms of the eigenenergies and eigenvectors of the 
energy matrix) depend on the following parameters: 
fermionic correlators C„ m (i — j) = (V'ra(i)V , m(j)) with 
|i — j| up to 4 lattice hops along the two main axis of 
the lattice and unknown bosonic correlators a M , 6 M , 
c M and d M (whose explicit expression can be found in 
Ref. [3]) and 
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a; = (ct"(i)c a (i)n(i)) - i (4^o - 1) (c^ (i)cr q c a (i)n q (i)) 
+ 4(25 M o-l)(c T (i)c i (i)c} a (i)c{' 3 (i)) (18) 

c ^ = ^ (c^^cfiOnfis)) - (4<5^o - 1) (c'W^cfiOri, (i B )) 
+ (2<5 M o - 1) (c T (i 5 )cj (j 5 )c{ (*i)c{(i)) 
+ \ (c\i)c(i 5 )n(h)} - — (45 M o - 1) ^c'^jffjcfe)!!, 
+ (25 M o - 1) (c T (ii)cj. (ii)c| (i 5 ) c |(i)) (19) 

d M = 4 (c'Mctiilnftis)) - — (4<5 M o - 1) (c f (»)<7,c (ii) n q (»i 3 )) 
+ (25 MO - 1) (c T (»is)cj (u 3 )c{ (ii)c|(i)) 

+ j (c t (i)c(ii 3 )n(ii)^ - ^ (4<5^o - 1) (J (i)a q c (ii 3 ) n q (ii)) 
+ (25„ - 1) (c T (ii) C| (ii) c\ (ii 3 ) cj(i)) (20) 

In order to compute the unknown bosonic correlators 
we can resort to the hydrodynamic constraints that re- 
quire that the bosonic spectra should be superlinear in 
momentum for long wavelength and that the suscepti- 
bility should be single valued at k = 0. Moreover, we 
can impose the local algebra constraint (n M (i)n M (i)) = 
(n) + 2(25 m0 - 1)A where D = (n}/2 - C 22 (0) is the 
double occupancy. According to this, we get four equa- 
tions which allow to compute the parameters a M , fe M , c M 
and d M in terms of a'^, and d'^. These latter param- 
eters, with respect to the no-primed ones, are made up 
of correlators containing operators centered at more 
distant sites. According to this, we expect that their 
values and relevance should be lower and in order to 
determine them we suggest the following decouplings 

< » 8 (2«5 m0 - 1) C a cc Ct - nC a cc - nC K cc (21) 

c m ~ 2 (2<5 M o — 1) C" c {Ccc + Cf c ) ~ 2^" c ~ 2^ c 

(22) 

df! ~ 2 (2<5 M o — 1) C cc (C cc + C' cc ) — ^Ccc — ~2^ cc 

(23) 

where k is the projector on the fourth nearest neigh- 
bors. 

In conclusion, we have reported a self-consistent 
scheme of calculations for the (spin and charge) bosonic 
sector of the t-t' Hubbard model. It is worth noticing 
that, within this scheme, the hydrodynamic constrains 
and the local algebra is preserved assuring that the 
known limits are conserved. Results of the presented 
scheme, easily attainable by solving numerically the 
self-consistent equations, will be presented elsewhere. 
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